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SECTION  I 


INTRODUCTION 

This  report  documents  the  third  version  of  a  computer  program,  the  Underwater  Shock  j 

i 

Analysis  (USA)  Code,  that  calculates  the  transient  response  of  a  totally  or  partially  sub-  I 

merged  structure  to  a  spherical  shock  wave  of  arbitrary  pressure-profile  and  source  loca- 

.  i 

tion.  The  structure  is  considered  to  be  linear-elastic  and  the  surrounding  fluid  is  treated  ' 
as  an  infinite  or  semi-infinite  acoustic  medium.  The  computational  model  for  the  structure 
is  constructed  through  the  use  of  an  auxiliary  discrete-element  (finite-element,  finite-  j 

difference)  code  of  choice  [1,2],  while  that  for  the  fluid  is  constructed  through  the  use  of 
the  Doubly  Asymptotic  Approximation  (DAA^)  or  the  Improved  Doubly  Asymptotic  Approximation  j 

(DAA2)  [3,4].  J 

As  implied  above,  this  manual  constitutes  a  revision  of  [5],  the  original  manual  for  the  ] 

USA  Code.  In  addition  to  various  maintenance  upgrades  introduced  into  the  code  since  it  was 
first  developed,  the  following  extensions  have  been  incorporated:  1)  Fluid  surface  elements 
for  wet-surface  segments  of  revolution  that  permit  a  circumferential  expansion  of  structural 
surface  motions  into  an  arbitrary  number  of  Fourier  components;  this  allows  the  inclusion  of 
both  beam-like  and  bar-like  motions  of  the  structure  as  a  special  case;  2)  An  imaging  formu¬ 
lation  that  extends  DAA  analysis  to  problems  involving  partially  submerged  structures  and 
structures  totally  submerged  near  the  fluid’s  free  surface;  3)  A  model  for  the  effects  of 
hulk  cavitation  on  incident-wave  excitation  that  provides  a  complete  description  of  fluid- 
particle  velocity  consistent  with  the  occurrence  of  surface  cutoff;  4)  The  introduction  of 
the  DAA2;  and  5)  Out-of-core  processing  of  the  fluid  equation  system. 

1.1  DOUBLY  SYMPTOTIC  APPROXIMATION 


The  principal  advantage  of  the  DAA^  and  DAA2  Is  that  they  model  the  acoustic  medium  sur¬ 
rounding  the  structure  as  a  membrane  covering  the  wet  surface  of  the  structure.  Hence  fluid 
motion  is  described  merely  in  terms  of  wet-surface  response  variables,  which  are  then  linked 
by  compatibility  relations  to  the  structural  response  variables.  Furthermore,  this  descrip¬ 
tion  is  a  simple  matrix  ordinary  differential  equation  with  desirable  compu tat ional  properties. 


The  principal  disadvantage  of  the  DAA^  is  that  it  constitutes  an  approximation  to  the 
"exact"  boundary-element  representation  of  the  surrounding  medium  [6,7].  The  DAA^  does  ap¬ 
proach  exactness  for  both  high-frequency  (early -time)  and  low  frequency  (late-time)  struc¬ 
tural  motions,  however,  and  effects  a  smooth  transition  between  the  two  asymptotes.  In  ad¬ 
dition,  it  has  exhibited  satisfactory  accuracy  in  a  variety  of  check  calculations  [4,6,8]. 
Hence,  in  view  of  its  desirable  computational  properties,  the  DAA^  is  considered  suitable 
for  engineering  analysis. 
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The  DAA^  is  an  improved  approximation  that  is  based  upon  the  DAA^;  however,  it  can 
describe  the  intermediate  frequency  range  better  than  the  DAA^ .  Computationally  it  has 
been  used  to  study  the  response  of  the  infinite  cylindrical  shell  under  a  plane  wave  step 
loading  in  which  significant  improvements  in  accuracy  have  been  noted.  Analytical  studies 
of  spherical  shell  response  [4]  also  show  such  marked  improvement. 

1.2  STAGGERED  SOLUTION  PROCEDURE 

The  governing  matrix  equation  for  structural  response  is  a  second-order  differential 
equation  in  time,  while  that  for  fluid  response  is  a  first-order  ordinary  differential  equa¬ 
tion.  Simultaneous  solution  of  these  equations  by  direct  step-by-step  numerical  integration, 
however,  is  unacceptably  expensive.  Hence  the  USA  code  utilizes  a  staggered  solution  pro¬ 
cedure  [9]  for  step-by-step  solution  of  the  equations  in  time. 

Now  a  staggered  solution  procedure  involves  a  response  extrapolation  at  each  time  step, 
which  usually  leads  to  numerical  instability  for  time  increments  exceeding  a  critical  value. 
Because  this  critical  value  may  be  unacceptably  small  for  many  computations,  the  governing 
equations  for  fluid  response  have  been  modified  in  such  a  way  that  unconditional  stability 
is  achieved.  Thus,  through  avoidance  of  both  direct  simultaneous  solution  and  conditional 
stability  constraints,  highly  efficient  computation  is  possible  for  the  greatest  variety  of 
cases. 

As  an  illustration  of  the  capabilities  of  USA,  a  transient  response  calculation  has 
been  performed  for  a  2490  degree-of-f reedom  (DOF)  structural  model  with  a  stiffness-matrix 
average  half-bandwidth  of  85  DOF.  The  central-processing-unit  (CPU)  time  on  a  Univac  1108 
required  for  the  280  time-step  calculation  (with  a  single  change  in  time  increment  during 
the  calculation)  was  28  minutes.  The  corresponding  time  on  a  CDC  6600,  on  which  the  code 
also  operates,  would  be  about  10  minutes. 

1.3  INPUT/OUTPUT 

The  USA  Code  requires  three  types  of  input  data  in  order  to  perform  its  function.  First, 
structural  mesh-geometry,  mass-matrix  and  stiffness-matrix  data  must  be  provided  by  the  struc¬ 
tural  analysis  code  used  by  the  analyst.  Second,  fluid  raesh-geometry  and  boundary-element 
data  must  be  furnished.  Finally,  location  of  the  fluid’s  free  surface,  charge  standoff,  inci¬ 
dent  pressure-profile,  and  time  integration  specifications  must  be  provided. 

The  code,  in  its  turn,  outputs  structural  displacement  and  velocity  histories,  and  fluid 
pressure  histories  for  the  wet  surface.  Response  data  post-processors  furnish  pseudo-velocity 
shock  spectra,  and  response-history  and  shock-spectrum  plots.  In  addition,  post-processors 
embedded  in  the  structural  analysis  code  may  be  used  to  obtain,  for  example,  stress  and  strain 
response  histories,  as  well  as  stress/strain-history  and  stress/strain-computer  plots.  As 
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currently  configured,  the  USA  Code  can  routinely  handle  problems  with  as  many  as 
3000  Structural  and  400  fluid  DOF  within  a  core  allocation  of  65000  decimal  words. 

1.4  SPECIAL  FEATURES 

A  number  of  special  features  are  incorporated  in  the  code.  First,  a  capability  has 
been  provided  to  handle  a  fluid  mesh  on  the  wet  surface  that  is  not  coincident  with  the 
surface  mesh  for  the  structural  model.  This  permits,  for  example,  the  use  of  a  refined 
structural  mesh  in  a  region  of  high  stress  gradients,  even  though  a  relatively  coarse  mesh 
is  retained  for  the  fluid. 

Second,  options  for  variable-increment  time  integration  and  computation  restart  are 
furnished.  The  former  allows  the  use  of  small  time  increments  during  periods  where  the  re¬ 
sponse  is  expected  to  be  varying  rapidly  in  time,  and  the  use  of  large  time  increments  for 
periods  characterized  by  a  slowly  varying  response.  The  latter  permits  the  division  of  a 
response  computation  into  segments,  so  that  the  analyst  may  examine  the  results  at  selected 
points  along  the  way.  Such  examination  is  facilitated  by  the  use  of  the  "printer-plot" 
routine  that  augments  the  usual  printout  data  with  response  plots  "drawn"  by  the  printer. 

Third,  the  code  incorporates  fluid  boundary  elements  for  both  general  and  body-of-re- 
volution  wet-surface  geometries  [10).  This  feature  is  especially  useful  for  compartment-by¬ 
compartment  analysis  of  a  submarine.  Such  an  analysis  utilizes  a  general-structure  discrete- 
element  model  of  a  particular  compartment  of  interest,  with  the  remainder  of  the  submarine 
modeled  as  a  bar/beam.  Hence  a  detailed  analysis  of  an  entire  submarine  may  be  performed 
with  several  discrete-element  models  of  moderate  size,  avoiding  the  use  of  a  single  gigantic 
model . 

Fourth,  the  analyst  can  use  either  the  DAA^  or  DAA^  through  input  of  a  single  scalar 
parameter  to  take  advantage  of  the  enhanced  accuracy  demonstrated  by  the  latter  in  computa¬ 
tions  for  idealized  geometries. 

Fifth,  out-of-core  processing  for  both  the  structural  and  fluid  equation  systems  frees 
the  user  from  concern  over  core  limitations  on  the  number  of  structural  and  fluid  elements 
in  his  model. 

Sixth,  free-field  shock  wave  input  to  the  structure  is  associated  with  a  spherical  wave 
and  can  be  input  for  use  with  linear  interpolation  or  cubic  spline  fitting  routines.  Pres¬ 
sure  histories  for  exponentially  decaying  incident  waves  can  be  automatically  generated. 

Fluid  pressure  and  particle  velocity  histories  cor responding  to  the  input  shock  are  displayed 
for  the  user  with  the  "printer-plot"  software  mentioned  above. 

Seventh,  routines  embedded  in  USA  can  be  used  to  facilitate  coupling  with  any  linear 
structural  analyzer  provided  lat  the  '  ffness  matrix  is  made  available  row  by  row  (or. 


J 


column  by  column).  These  routines  assemble  the  matrix  in  the  partitioned  skylined  format 
required  by  USA.  In  addition,  the  stiffness  matrix  can  be  displayed  for  check  purposes. 

1.5  SPECIAL  NOW- FEATURES 

Some  features  of  modest  complexity  have  yet  to  be  incorporated  into  the  USA  Code. 
First,  an  option  for  automatic  time-step  integration  would  free  the  analyst  from  having  to 
select  integration  time  increments  in  accordance  with  his  expectations  regarding  response 
behavior.  Second,  the  ability  to  handle  very  large  problems  would  be  useful  in  those  cases 
where  structural  segmentation  is  not  possible.  Third,  a  capability  to  treat  banded  struc¬ 
tural  mass  and  damping  matrices  would  be  helpful,  in  order  to  accommodate  structural  ana¬ 
lyzers  that  produce  such  matrices.  Fourth,  a  means  to  handle  the  matrices  produced  by 
acoustic  elements  based  on  a  pressure  formulation  [11,12]  is  desirable,  in  that  such  ele¬ 
ments  permit  highly  efficient  modeling  of  internal  fluid  volumes.  Finally,  a  capability  to 
treat  localized  nonlinear  structural  behavior,  such  as  that  exhibited  by  nonlinear  equip¬ 
ment  mounts,  would  permit  highly  efficient  treatment  of  such  behavior;  for  widespread  non¬ 
linear  effects,  however,  recourse  must  be  made  to  a  fully  nonlinear  code,  such  as  USA- 
STAGS  [13]. 

An  important  feature  of  greater  complexity  has  yet  to  be  introduced  into  the  code. 

This  is  a  treatment  of  hull  cavitation,  which  may  substantially  affect  structural  response 
for  incident  shock  waves  of  short  duration.  The  introduction  of  this  feature  requires  the 
accurate  treatment  of  highly  nonlinear  phenomena,  and  presents  a  challenging  task  for 
future  work. 


SECTION  II 


THEORY 

This  section  describes  the  theoretical  foundation  of  the  USA  Code.  It  is  constructed 
as  an  overview,  with  coverage  of  details  left  to  referenced  papers  and  reports. 


2.1  STRUCTURAL  RESPONSE  EQUATION 

The  matrix  ordinary  differential  equation  for  the  dynamic  response  of  a  linear-elastic 
structure  is  f 1 | 


Mx  +  Cx  +  Kx=f  (2.1) 

~s  —  ~s  —  ~s  —  — 

where  x  is  the  structural  displacement  vector,  ,  Cg  and  are  the  structural  mass, 

damping  and  stiffness  matrices,  respectively,  is  the  external  force  vector,  and  a  dot 

denotes  a  temporal  derivative.  Generally,  and  are  highly  banded,  symmetric 

matrices  of  large  order;  at  present,  the  USA  Code  considers  to  be  diagonal  and  C  to 

be  zero. 

For  excitation  of  a  submerged  structure  by  an  acoustic  wave,  .f  is  given  by 

1  =  -£Af  (£].  +  £s)  (2.2) 

where  and  are  nodal  pressure  vectors  for  the  wet-surface  fluid  mesh  pertaining  to 

the  (known)  incident  wave  and  the  (unknown)  scattered  wave,  respectively,  A^  is  the  diag¬ 
onal  area  matrix  associated  with  elements  in  the  fluid  mesh,  and  G  is  the  transformation 
matrix  that  relates  the  structural  and  fluid  nodal  surface  forces.  More  will  be  said  about 
G  in  the  next  subsection. 

2.2  DAA1  EQUATION 

The  Doubly  Asymptotic  Approximation  may  be  written  [3,4] 

Mf  is  +  PcAf£s  =  PcMfus  (2.3) 

where  u^  is  the  vector  of  scattered-wave  fluid-particle  velocities  normal  to  the 
structured  wet  surface.  p  and  c  are  the  density  and  sound  velocity  of  the  fluid, 
respectively,  and  is  the  symmetric  fluid  mass  matrix  for  the  wet-surface  fluid  mesh. 

This  matrix  is  produced  by  a  boundary-element  treatment  of  Laplace1 S  equation  for  the 
irrotational  flow  generated  in  an  infinite,  inviscid,  incompressible  fluid 
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by  motions  of  the  structure's  wet  surface;  it  is  fully  populated  with  non-zero  matrix  ele¬ 
ments.  When  transformed  into  structural  coordinates,  the  fluid  mass  matrix  yields  the  added 
mass  matrix,  which,  when  combined  with  the  structural  mass  matrix,  yields  the  virtual  mass 
matrix  for  motions  of  a  structure  submerged  in  an  incompressible  fluid  [14]. 

As  mentioned  in  Section  I,  the  approximate  relation  (2. 3)  is  called  "doubly  asymptotic" 
because  it  approaches  exactness  in  both  the  high-frequency  (early-t iine)  and  low-frequency 
(late-time)  limits.  For  high-frequency  motions,  |£g|  >>  | pg |  ,  so  that  (2 . 3)  approaches  the 

relation  £g  =  pcu^  ,  which  is  the  correct  limit  for  short  acoustic  wavelengths.  For  low- 
frequency  motions,  |  <<  |pQ |  ,  so  that  (2. 3) approaches  the  incompressible-flow  relation 

jig  ,  which  is  the  correct  limit  for  long  acoustic  wavelengths. 

For  excitation  by  an  incident  acoustic  wave,  is  related  to  structural  response  by 

the  kinematic  compatibility  relation 

GT  x  -  u. !  +  u g  (2.4) 

where  the  superscript  "T"  denotes  matrix  transposition.  Equation  £.4}  expresses  the  con¬ 
straint  that  normal  fluid-particle  velocity  match  normal  structural  velocity  on  the  wet 
surface  of  the  structure.  The  fact  that  the  transformation  matrix  relating  those  velo- 
cities  is  G  follows  from  the  invariance  of  virtual  work  with  respect  to  either  of  the 
wet  surface  coordinate  systems.  Generally,  G  is  a  rectangular  matrix  whose  height 
greatly  exceeds  its  width,  inasmuch  as  the  number  of  structural  DOF  usually  exceeds  con¬ 
siderably  the  number  of  fluid  DOF. 

2.3  INTERACTION  EQUATIONS 

The  introduction  of  (2.2)  into  (2.1)  and  (2.4)  into  (2.3)  yields  the  interaction  equa¬ 
tions 

M  x  +  C  x  +  K  x  -  -  G  A-  (pT  +  Pc) 

-  S  —  ~S  -  ~3  —  ~  1 

T  (2.5) 

Mf  Ps  +  P  c  Af  ps  =  p  c  Mf  (G  x  -  uL) 

These  equations  may  be  solved  simultaneously  at  each  time  step  by  the  transfer  of  -  G 
and  p  c  M  G  x  to  the  left  sides  of  their  respective  equations.  Such  a  procedure  is 
exceedingly  expensive,  however,  because  of  the  large  connectivity  of  the  coefficient  matrix 
involved.  As  mentioned  in  Section  I,  efficient  computation  is  possible  through  the  appli¬ 
cation  of  a  staggered  solution  procedure  that  is  unconditionally  stable  with  respect  to 
the  choice  of  time  increment. 

The  simplest  implementation  of  the  staggered  solution  procedure  recommended  in  [9]  may 
be  effected  as  follows.  M  is  taken  to  be  diagonal  and,  to  allow  for  the  possibility  that 
M  may  have  zero  entries  for  rotational  DOF,  0  is  constructed  such  that  onlv  the  transla- 
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tional  DOF  for  the  structure  couple  with  the  fluid  DOF  [see  (2.4)];  then  the  first  of  (2.5) 

T-- 

may  be  partitioned  to  obtain  G  x>  which  may  then  be  introduced  into  the  second  of  (2.5). 
Premultiplication  of  the  resulting  equation  by  1  .  „  1 


AF  Mr 
pc  ~f 


then  yields 


J. 

pc 


-  p. 


s  +  <8n+  Ba>  £S 


-if  5T  &l(J,  i  +  k,  x) 


-<Ss  El  +  if  "l> 


(2.6) 


-1  X  -1  -1 

where  D^  =  and  Dg  =  A^  G  Mg  G  are  symmetric ,  and  where  Mg  is  a  diagonal 

matrix  with  each  nonzero  element  given  as  the  reciprocal  of  the  corresponding  nonzero  element 

of  Mg  and  each  zero  element  mirroring  the  corresponding  zero  element  of  Mg  .  The  first 

of  (2.5)  and  (2.6)  are  herein  termed  "the  augmented  interaction  equations". 


2.4  SPHERICAL  INCIDENT  WAVE 


Each  element  of  the  vectors  and  for  a  spherical  incident  wave  are  given  by 


/  -s  \ 

-  r7  pi  (c  *  -V-) 

“ii/o  =  l  £  Putt)  +  ^7  y1 


(2.7) 


where  S  is  the  "charge  standoff",  i.e.,  the  distance  between  the  origin  of  the  incident 
spherical  wave  and  the  nearest  point  on  the  structured  wet  surface,  is  the  distance 

from  the  origin  of  the  incident  spherical  wave  to  the  ith  fluid  node  on  the  wet  surface, 
is  the  cosine  of  the  angle  between  the  vector  corresponding  to  and  the  wet-surface 

normal  at  the  ith  fluid  node,  and  p^-(t)  is  the  incident-wave  pressure-profile  defined  at 
=  S  .  For  a  shock  wave,  p^(t)  is  discontinuous  at  t  =  0  and  the  u^(t)  contain 
singularities . 

In  order  to  remove  shock-wave  singularities  from  u^  in  (2.6),  a  modified  pressure 
vector  is  defined  as 


%  =  ££l  +  £s  (2.8) 

where  _T  is  a  diagonal  matrix  with  direction-cosine  elements  .  The  introduction  of 

(2.8)  into  (2.6)  and  the  first  of  (2.5),  followed  by  utilization  of  the  second  of  (2.7)  then 
vie  Ids  tile  modified,  augmented,  interaction  equations 


M 


x  +  Cg  x  + 


x 


J__ 

cc 


A  Em  +  (2fl+  Ss  }  Em 


-  5  ~f  [Em  +  (-r  -  £i] 

'  A f  £T  iQ1  (C„  x  +  Ks  x)  -  _Hp L 


(2.9) 
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in  which  I  is  the  identity  matrix,  and 


D  -  (D  +  D  -  —  Ac  R  X)  T 
~s  1  p  ~  ~ 


(2.10) 


where  -R  is  the  diagonal  matrix  formed  by  the  distances  R^.  Equations  (2.9)  (with 
jC^  -  0)  are  the  equations  solved  by  the  USA  Code  to  determine  the  structural  responses 
x  and  x,  and  the  wet-surface  pressures  £  =  (1^  -  V)  4- 

2.5  FREE  SURFACE  EFFECTS 

When  a  structure  is  partially  submerged,  or  when  a  totally  submerged  structure  lies 
near  the  free  surface  of  a  semi- inf inite  fluid,  imaging  techniques  may  be  utilized  to  en¬ 
sure  that  the  total  pressure  vanishes  at  the  free  surface.  (This  implies  that  the  effects 
of  gravity  are  negligible  in  this  class  of  problems,  which  they  generally  are.)  In  this 
case,  the  interactive  system  consists  of  an  infinite  fluid  domain,  the  structure  S+,  and 
its  image  S_  (see  Figure  2-1).  The  incident  wave  now  consists  of  a  (positive)  primary 
wave  plus  a  (negative)  image  wave,  the  latter  emanating  from  the  image  of  the  primary 
wave’s  origin.  Zero  pressure  at  the  free  surface  is  therefore  maintained  if  the  motions 
of  S_  are  constrained  to  be  opposite  to  those  of  S+. 

The  kinetic  energy  T  ,  the  Rayleigh  dissipation  function  D  ,  the  potential  energy  V  , 

s  s  s 

and  the  work  potential  ng,  for  the  structural  system  S+  +  S_  are  given  by 

T  "  i  (*+  JtL  5+  +  JL  *-> 


i  Cg  X+  +  X_  Cs  x.) 


i  <£  K  5+  +  xl  K  x_) 


(2.11) 


T  T 

ns  -  -x+  f+  -  X-  f- 

The  appropriate  constraints  are  x_  = -x+  and  f_  =  -f+,  so  that  (2.11)  become 

T  -  xT  M  icf 

D  =  x!  C  x± 

s  —  i  ^ 


(2.12) 


V  =  x+  K  x+ 
s  ”  — r 


nQ  =  -2  1+ 


The  DAA  kinetic  energy  and  work  potential  for  the  fluid  system  may  be  written 


where 


Tf  ^  Hs  ~f  Hs 

*T  1  *t 

nf  =  BS  - s  ~f  E< 


M  M' 


(2.13) 


(2.14) 


and  an  asterisk  denotes  temporal  integration.  The  submatrix  M  accounts  for  added  mass 
coupling  between  wet-surface  elements  on  S+  and,  similarly,  between  elements  on  S_;  Mf 
accounts  for  added  mass  coupling  between  elements  on  S+  and  elements  on  S"_.  The  con¬ 
straints  for  the  fluid  system  are 


V  *  2s 


(2.15) 


so  that  (2.13)  becomes 


Tf  =  Hs+ 


(2.16) 


-2  &  [A  ps+  +  i  (M  -  M-)  ps+] 


The  application  of  Lagrange’s  equation  [15]  to  (2.12)  and  (2.16)  now  yields 

M  x  +  C  x  +  K  x,  =  f 
~s  -+  -+  s  -+  -+ 


—  (M  -  M')  p  +  A  p. 


(2.17) 


Qi  -  S')  hs+ 


Also,  (2.2)  and  (2.4)  must  be  modified  to  include  the  effects  of  both  the  incident  prim¬ 
ary  and  image  waves.  This  gives 


1+  =  -£  A  (El+  +  £I+  +  Es+) 


~  *+  =  +  Bx+  +  BS+ 


(2.18) 


where,  e.g. ,  denotes  incident-wave  pressure  on  S+  associated  with  the  image  wave. 
The  introduction  of  (2.18)  into  (2.17)  then  yields  the  doubly  asymptotic  interaction 
equations  for  problems  involving  a  free  surface 
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(2.19) 


M  x :  +  C  x .  +  K  x  =  G  A  (p*  +  pT,  +  p„  ,  ) 

~s  -+  ~s  -+  ~  ~  M+  -1+  -S  + 

(51  *  51)  Es+  +  Pc  A  £s+  =  pc  (M  -  MT)  (GT  x+  -  u++  -  u*+) 

A  comparison  of  (2.19)  with  (2.3)  reveals  that  the  effects  of  the  free  surface  are  embod¬ 
ied  in  the  image-wave  pressure  and  fluid-particle-velocity  vectors,  and  in  the  modified 
added-mass  matrix. 

Finally,  augmentation  of  (2.19)  to  secure  unconditional  stability,  followed  by 
introduction  of  the  modified  pressure  [cf.  (2.8)] 

-M  =  ~+  2l+  +  ~+  2l+  +  Bs+  (2.20) 

to  remove  shock-wave  singularities,  proceeds  as  described  in  Subsections  2.3  and  2.4. 

The  modified,  augmented  interaction  equations  corresponding  to  (2.9)  for  the  infinite 
f  luid  medium  are  then  readily  obtained. 

It  is  important  to  mention  at  this  point  that  the  DAA  formulation  just  described 
does  not  account  for  high-frequency  scattered  waves  from  S_  that  impinge  upon  S+.  For 
most  floating  structures,  such  waves  are  not  generated,  as  the  wet  surfaces  of  and  S 
usually  intersect  to  form  a  convex  surface;  they  are  generated,  however,  for  a  totally 
submerged  structure  lying  near  the  free  surface.  Even  so,  it  has  been  shown  that,  as 
far  as  structural  response  is  concerned,  the  effects  of  the  scattered  wave  are  generally 
negligible  [16].  In  other  words,  the  response  is  basically  driven  by  the  incident 
primary  and  image  waves. 

2.6  BULK  CAVITATION 

In  the  absence  of  bulk  cavitation,  the  imaging  method  serves  as  a  useful  device  to 
model  the  reflection  of  free-field  waves  from  the  fluid's  free  surface.  The  occurrence 
of  bulk  cavitation  near  the  surface,  however,  changes  that  simple  acoustic  reflection 
problem  into  a  complex  reflection-refraction  problem,  as  indicated  in  Figure  2-2.  If 
refraction  distortions  produced  by  a  relatively  thin  cavitated  region  are  not  too  severe, 
however,  bulk  cavitation  effects  will  still  appear  to  the  structure  as  emanating  from  an 
image  source. 

Experimental  records  of  free-field  pressure  histories  for  compact  charges  exhibit 
the  behavior  shown  in  Figure  2-3  [17].  The  dashed  line  denotes  the  history  produced  by 
a  negative- image  model,  while  the  horizontal  line  indicates  that  the  effect  of  bulk  cavi¬ 
tation  is  to  "cut  off"  the  pressure  at  a  cavitation  threshold.  The  approximate  treatment 
introduced  here  involves  p re- examination  of  the  image-based  free-field  pressure  at  the 
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Figure  2-3  Free-Field  Pressure  as  a  Result 
of  Free-Surface  Reflection 
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Figure  2-4  Geometry  of  Primary  and  Image  Waves 
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standoff  point,  i.e.,  the  point  on  the  submerged  structure  closest  to  the  charge.  When¬ 
ever  that  pressure  becomes  negative  to  the  extent  that  its  magnitude  exceeds  the  absolute 
ambient  pressure  at  that  depth,  a  positive  contribution  is  incorporated  into  the  negative- 
image  source  so  that  the  free-field  pressure  at  the  standoff  point  never  dips  below  the 
cavitation  threshold.  When  the  structure’s  overall  dimensions  are  small  relative  to  the 
distance  from  the  structure  to  the  cavitated  region,  the  effects  of  the  positive  contri¬ 
bution  will  not  vary  appreciably  in  the  vicinity  of  the  structure. 


The  preceding  discussion  leads  to  the  following  development.  The  two-source  model  of 
Figure  2-4  yields  as  the  free-field  pressure  at  any  point  P 


pp(t) 


(t 


R .  -  S  ,,  R  -  S 

-V->  +  rp- (t  '"V-> 


(2.21) 


where  p+(t)  =  p  (t)  =  0  for  t  <  0.  At  the  standoff  point,  (2.21)  becomes 

s  -  s 

pQ(t)  -  p+(t)  +  -—  p_  (t  -  ~ —  ) 


(2.22) 


Now  p  (t)  =  -p+(t)  as  long  as  the  resulting  p^(t)  exceeds  the  cavitation  threshold  so  that 
"cutoff"  does  not  occur;  otherwise,  p^(t)  remains  at  the  threshold  value  ~(PA  +  Y  Z) 
where  p^  is  atmospheric  pressure  and  y  is  the  fluid's  weight  density.  Hence,  during  the 
"cutoff  period", 

s  s  -  s 

P_(t)  -  -  -f  fp+  <t  +  - - )  +  PaHZ]  (2.23) 

The  model  just  described  fits  a  prescribed  free-field  pressure  history  at  the  stand¬ 
off  point  in  such  a  way  that  surface  cutoff  effects  appear  to  the  structure  as  emanating 
from  an  image  source.  Because  the  model  is  complete,  it  also  provides  the  free-field 
fluid-particle-velocity  information  required  for  DAA  calculations.  The  usefulness 
of  the  model  has  been  demonstrated  from  the  results  of  free-field  tests  specifically 
designed  to  produce  both  pressure  and  fluid-particle-velocity  data  I 18 j . 

2 . 7  usa-daa2  IMPLEMENTATION 

The  Improved  Doubly  Asymptotic  Approximation  [DAA^]  can  be  written  as  [4] 


Mf£s  +  pc6f£s  +  Pc3f  Af£g  = 
pc  llJ  f(QT'x- u  j)  £f  Mf(GTx-u  r)J 
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where 


Qf  =  HpcAf  Mj1  (2.25) 

All  vector  and  matrix  quantities  in  the  above  are  related  to  the  same  finite  element 
wet-surface  fluid  mesh  as  that  used  for  the  lowest  order  DAA  (DAA^)  described  in 
Section  2.2. 

Note  that  DAA^  is  a  second-order  equation,  whereas  DAA^  is  a  first-order  equation.  In 
addition,  DAA^  includes  a  new  scalar  parameter  n  that  appears  in  (2.25).  It  can  be 
established  from  physical  considerations  1 that  must  be  bounded  as 


0  <  <  i 


(2.26) 


A  precise  choice  of  n  is  apparently  mt  pLedicated  by  any  fundamental  principle. 
Hence  it  must  be  regarded  at  this  time  as  a  factor  which  may  be  adjusted  to  achieve 
optimum  accuracy  for  a  particular  problem.  In  [4],  it  is  observed  that  T)  »  1  leads 
to  the  best  accuracy  for  a  spherical  shell. 


In  order  to  implement  DAA~  (2.24)  is  first  integrated  once  in  time  and  multiplied 
-1  1 

through  by  .  Equation  (2.25)  is  then  substituted  into  the  result  and  a  new 

variable,  the  scattered  pressure  integral  ,  is  defined  by 


* 

£  < 


(2.27) 


where  an  asterisk  denotes  temporal  integration.  The  result  is 


2  2 


£f  Ss  +  pc5fiis+  np  c  5 f 2  as  = 


T  2  2  t 

pc  A  ^  (G  x  -u  1)  +  r)  P  c  D  (G  x-u_j.) 


(2.28) 


where 


D 


f  2 


(2.29) 
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It  is  noted  that  (2.28)  is  symmetric  and  that  D  ^  has  already  been  defined  following 

(2.6). 

To  avoid  shock-wave  singularities  in  u  ,  the  relation  for  a  spherical  shock  is 
used  as 


pcuI  =  r  (£  j  +  c  R  £];) 


(2.30) 


while  the  modified  pressure-integral  vector  is  defined  as 


aM  =  + 


(2.31) 


Substitution  of  (2.30,  (2.31)  into  (2.28)  then  gives 


AfqM+Pc2fl^M  +  nP  C  5f2^M 

PC  Af  GT  X  +  n  P2  C2  Df]GT  X  +  c  [(1-n)  P  Dfl  -  A  f  R_1]  [£l 


(2.32) 


9  2  1  -1  * 

+  n  pc  (Df2  -  ±  Dfl  R  )  r  Px 


where  the  identity 


r'1  r  =  r  r  1 


(2.33) 


has  been  used  in  (2.32),  as  both  matrices  are  diagonal.  Associated  with  (2.32)  is 
the  structural  equation  of  motion 


M  k  +  C  x  +  K  x  *  -C  A-  [q„  +  (I -T  )  p  T] 
s  —  ~  s  ~  s  ~  ~  i  M  ■*- 1 


(2.34) 


Equations  (2.32)  and  (2.34)  define  the  DAA^-mod if ied  interaction  equations  that  are 
solved  according  to  the  staggered  solution  strategy;  hence  an  examination  of  stability 
must  be  conducted.  It  has  been  shown  that  the  step-by-step  integration  of  (2.32)  and 
(2.34)  is  conditionally  stable;  however,  no  systematic  study  of  stability  has  yet  been 
undertaken.  In  view  of  the  fact  that  unconditional  stability  was  achieved  for  USA- 
DAA^  by  augmentation,  and  that  (2.24)  is  essentially  the  DAA^  with  a  correction  term, 
augmentation  of  (2.32)  was  carried  out  in  the  same  manner  as  that  used  for  DAA^ . 
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Accordingly,  (2.34)  is  first  solved  for  x  and  substituted  into  (2.32)  to  give 

!fa»*  "  (2  n+  2  8>  Am  +  n  p2  =2  5f2  a.  " 

-  PcAfGTM-sV,i  +  Ka*>  +np2c2DflGT  i 

(2.35) 

-  pc  {?s-  [ds+  d-n)Dfl  -  i  Af  R_1]  r}  £j 
+  n  p2  c2  (d^Jd^r-1)  T  £i 

where  Dq  has  already  been  defined  following  (2.6). 

Equations  (2.34)  and  (2.35)  are  the  DAA2“modif ied,  augmented  interaction  equations 
that  have  been  implemented  in  the  USA  Code. 
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SECTION  III 


ORGANIZATION 

The  USA  Code  has  been  written  in  standard  FORTRAN  IV  for  use  on  both  Univac  and  CDC 
computers.  Machine  dependency  has  been  isolated  in  one  utility  program  described  below. 
Program  modularity  has  been  strictly  enforced,  with  communication  between  computational 
modules  controlled  by  means  of  a  data  management  system. 

The  basic  structure  of  the  code  is  shown  in  Fig.  3-1.  The  structural  preprocessor  is 

a  separate  code  selected  by  the  user  to  provide  the  computational  model  for  the  structure. 

The  skyline  utility  merely  reformats  M  and  K  as  provided  by  the  structural  preprocessor 

for  processing  by  the  USA  Code  (recall  that  C  is  taken  as  zero) .  The  fluid  mass  prepro- 

~  s 

essor  forms  A^,  M^,  D^,  D^  and  G  using  a  virtual  memory  simulator  for  out-of-core  process- 

~  ~  ~  T  -  \ 

ing,  while  the  matrix  augmentation  preprocessor  forms  D  and  AG  M  [see  (2.9) J.  The 

~  s  i  ~  >v  s 

main  processor  is  the  time  integrator,  which  forms  T  and  H  and  then  solves  (2.9)  in  step- 
by-step  fashion  using  the  staggered  solution  procedure.  The  response  postprocessor  pro¬ 
vides  tabular  and  graphic  output  for  the  computed  kinematic  responses  as  well  as  pseudo¬ 
velocity  shock  spectra.  Finally,  the  data  manager  controls  the  flow  of  data  between  pro¬ 
cessors.  More  detailed  descriptions  of  the  various  program  components  follow,  while  infor¬ 
mation  required  for  utilization  of  the  code  is  contained  in  Appendices  A  through  D. 

3.1  TEE  DATA  MANAGER  DMGASP 

DMGASP  is  a  self-contained  utility  module  that  functions  as  a  manager  of  auxiliary  stor¬ 
age  and  as  the  focal  point  for  all  block  input/output  activities  [19].  Constituting  the 
lowest  level  of  the  NOSTRA  Data  Management  System  [20  ] ,  it  carries  out  the  direct  transfer  of 
data  blocks  between  core  and  peripheral  storage.  (The  terminology  "direct  transfer"  is  used 
here  to  denote  unformatted  and  unbuffered  data  transmission.)  The  basic  auxiliary  storage 
management  operations  embodied  in  DMGASP  are 

•  Activate  storage  device 

•  Position  device 

•  Read  data  block  from  device 

m  Write  data  block  on  device 

•  Deactivate  device 

In  the  USA  Code,  DMGASP  is  operated  as  a  stand-alone  I/O  package  that  receives  direc¬ 
tives  directly  from  the  master  processors.  Assembly  language  versions  of  DMGASP  currently 
exist  for  UNIVAC  1100  EXEC-8,  CDC  SCOPE  3.4  (NOS/BE),  and  CDC  NOS  operating  systems;  hence 
the  USA  Code  may  be  used  only  on  these  systems  at  this  time. 
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3.2  THE  VIRTUAL  MEMORY  SIMULATOR  VMSYST 


VMSYST  is  a  virtual  storage  simulator  for  computers  that  are  not  built  around  a  vir¬ 
tual  memory  system  [21].  All  data  in  the  virtual  system  is  partitioned  into  pages,  which 
are  blocks  of  consecutive  data  words  of  a  fixed  page  size.  Pages  residing  in  core  in  the 
page  buffer  are  called  active  pages.  Inactive  pages  are  resident  in  auxiliary  storage 
only.  In  this  utility  the  page  and  page  buffer  sizes  can  be  conveniently  adjusted  to  suit 
the  application.  Input  and  output  to  auxiliary  storage  is  handled  by  DMGASP;  otherwise 
VMSYST  is  written  in  transportable  FORTRAN. 

The  primary  advantage  of  a  virtual  memory  system  is  the  efficient  processing  of  many 
small  records  such  as  columns  or  rows  of  large  full  matrices  that  can  be  treated  as  vec¬ 
tors.  In  essence  VMSYST  keeps  track  of  whether  a  desired  block  of  data  is  resident  in 
core  in  the  page  buffer,  or,  has  been  moved  to  an  external  storage  device  by  DMGASP.  If 
it  is  not  currently  resident  in  the  page  buff er,  VMSYST  retrieves  it  and  makes  it  available 
to  the  application  program.  This  double  movement  of  data  is  the  major  price  paid  for  the 
benefits  of  the  virtual  system. 

In  USA,  VMSYST  is  used  for  the  out-of-core  generatiw  *  of  the  fluid  mass  matrix  des¬ 
cribed  in  Section  3.5. 

3.3  THE  STRUCTURAL  PREPROCESSOR 

This  is  a  user-provided  code  that  assembles  the  structural  mass  and  stiffness  matrices 
and  generates  information  that  relates  the  internal  and  external  descriptions  of  the  struc¬ 
tural  DOF.  Input  typically  includes 

•  Mesh  geometry 

Coordinate  systems 

Node  locations 

•  Element  definitions 

Type 

Connectivity 

•  Material  properties 

Mass  density 

Moduli 

•  Constraints 

Symmetry  conditions 

Element  external  constraints 

Element  internal  constraints 
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Fluid  internal  to  the  submerged  structure  must  be  included  in  the  structural  model. 

At  this  time,  USA  treats  only  diagonal  mass  matrices  associated  with  a  lumped  mass  repre¬ 
sentation  of  the  structure,  and  only  single  precision  matrices  can  be  processed. 

3.4  THE  SKYLINE  UTILITY 

This  preprocessor  converts  the  structural  mass  and  stiffness  matrices  generated  by 
thi  structural  preprocessor  into  the  internal  "skyline"  format  required  by  the  USA  time 
integration  processor  [22,23].  As  there  are  a  variety  of  ways  to  store  large,  sparse, 
symmetric  matrices,  virtually  any  structural  preprocessor  that  is  to  be  used  with  the  USA 
Code  will  require  a  utility  package  to  change  the  storage  format.  At  this  time,  conver¬ 
sion  utilities  have  been  written  for  SPAR  [24],  NASTRAN  [25]  and  GENSAM  [26].  User  in¬ 
structions  for  constructing  the  skyline  utility  for  other  structural  codes  are  given  in 
Appendix  E.  As  is  noted  there,  USA  now  contains  a  set  of  utilities  to  facilitate  this 
process . 

Figure  3.1  shows  2  paths  to  mass  storage  from  the  skyline  utility.  The  SPAR  converter 
uses  DMGASP  for  both  input  and  output,  whereas  the  NASTRAN  converter  uses  unformatted 
buffered  FORTRAN  commands  for  input  and  DMGASP  for  output. 

Constraints  are  also  handled  differently  in  these  two  utilities.  NASTRAN  provides 
a  reduced  stiffness  matrix  which  already  incorporates  any  prescribed  constraints.  SPAR 
does  not;  however,  USA  has  the  ability  to  apply  constraints  due  to  symmetry  or  attachment 
to  ground  during  the  time  integration.  Structural  DOF  that  must  be  set  to  zero  ate 
flagged  by  the  skyline  utility  [23]. 

3.5  THE  FLUID  MASS  PREPROCESSOR  FLUMAS 

This  code  constructs  the  fluid  mass  matrix  for  a  structure  submerged  in  an  infinite, 
inviecid,  incompressible  fluid  by  the  boundary  element  technique  [14],  In  addition, 

FLUMAS  can  form  the  mass  matrix  for  a  body  in  the  vicinity  of  a  fluid  free  surface  through 
the  use  of  imaging  techniques.  Out-of-core  processing  is  facilitated  by  use  of  the  vir¬ 
tual  memory  system  VMSYST  so  that  core  size  is  not  a  limitation  on  the  number  of  fluid  DOF. 
The  code  also  generates  fluid  mesh  data  and  a  set  of  transformation  coefficients  relating 
the  structural  and  fluid  DOF.  The  computation  of  these  coefficients  is  based  upon  the 
use  of  centroidal  nodes  for  the  fluid  elements  and  the  assumption  of  a  bilinear  variation 
of  displacement  over  the  surface  of  each  structural  element.  This  assures  that  the  des¬ 
cription  of  the  fluid  pressure  forces  in  the  two  mesh  systems  is  statically  equivalent 
without  inducing  moments  at  the  structural  nodes.  Finally,  the  code  generates  the  sym¬ 
metric  matrices  and  g^  t*vat  appear  in  the  computational  form  of  the  DAA^  and  DA^ 

equations  which  involve  the  inverse  of  the  fluid  mass  matrix. 
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FLUMAS  contains  a  refined  formulation  for  the  fluid  mass  matrix  that  includes  the 
primary  effects  of  element  curvature.  In  addition,  it  has  the  capability  to  treat  struc¬ 
tures  containing  both  general  geometry  and  arbitrary  axis,  multi-branch,  multi-harmonic 
surface-of-revolution  components,  as  described  in  [10].  The  code  can  also  efficiently 
construct  the  fluid  mass  matrix  for  a  body  with  one  or  two  planes  of  symmetry  by  using  a 
mesh  which  covers  1/2  or  1/4  of  the  surface.  Symmetric  or  anti-symmetric  fluid  motions 
can  then  be  imposed  on  the  portions  of  the  surface  not  covered  by  the  mesh.  Two- 
dimensional  "plane-strain"  behavior  of  long  cylinders  can  also  be  simulated.  The  code 
contains  an  automatic  mesh  generator  for  cylindrical  surfaces  and  an  improved  error  exit 
control  that  facilitates  fluid  mesh  debugging.  Finally,  a  useful  diagnostic  tool  in  the 
code  is  the  capability  to  solve  the  fluid-boundary-mode  problem  M^_u  =  XA^ii  [14]. 

Typical  input  data  for  this  processor  includes 

•  Mesh  geometry 

Fluid  Wet-Surface  Mesh 
Structure  Wet-Surface  Mesh 

•  Element  definitions 

General  curved  surface 
Surface  of  revolution 

•  Material  property 

Mass  density 
DAA2  parameter 

•  Constraints 

Location  of  free  surface 
Half  model 
Quarter  model 
Long  cylinder 

Node  reassignment  in  fluid-structure  transformation 
A  detailed  description  of  the  required  input  data  is  given  in  Appendix  A. 

3.6  THE  AUGMENTED  MATRIX  PREPROCESSOR  AUGMAT 

This  processor  accepts  data  from  the  structural  and  fluid  analyzers  to  construct  the 
specific  matrices  required  for  solution  of  the  augmented  Eqs.  (2.9)  or  (2 .34) / (2 . 35) . 

The  output  of  this  code  includes  not  only  the  required  matrices  in  skyline  form,  but  also 
a  distillation  of  the  output  from  both  the  structural  and  fluid  processors.  This  has 
been  done  so  that  only  one  permanent  file  need  be  referenced  as  input  to  the  time  inte¬ 
grator;  this  results  in  improved  data  handling  and  core  usage.  In  contrast  to  earlier 
versions  of  USA,  AUGMAT  does  not  form  the  fluid  matrices  D  ^  and  D  f ^  but  rather  puts 
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them  in  the  skyline  format  required  by  SKYPUL.  Dfl  and  D  are  now  formed  only  in 
FLUMAS.  The  structural  stiffness  matrix  can  also  be  displayed  in  AUGMAT  for  checkout 
purposes.  Input  to  this  code  involves  the  following  information 

9  Mass  matrices 
Fluid 
Structure 

•  Structural  DOF  correspondence  table 

External  and  internal  node  descriptions 

Factorization  order 

DOF  reduction  due  to  constraints 

•  Fluid  mesh  geometry 

Global  coordinates  of  fluid  nodes 
Direction  cosines  for  nodal  surface  normals 
Areas  of  fluid  elements 

9  Fluid/structure  DOF  transformation  coefficients 

9  Fluid  material  properties 
DAA^  parameter 

•  Constraints 

Half  model 
Quarter  model 

Although  this  constitutes  a  substantial  amount  of  information,  almost  all  of  it  is  re¬ 
trieved  from  permanent  data  files.  A  detailed  discussion  of  the  required  input  data  is 
contained  in  Appendix  B. 

3.7  THE  TIME  INTEGRATION  PROCESSOR  TIMINT 

This  main  processor  constitutes  an  implementation  of  the  unconditionally  stable 
staggered  solution  technique  developed  in  [9]  for  DAA^.  The  primary  output  is  a  set  of 
permanent  data  files  that  contain  nodal  values  for  structural  displacement,  structural 
velocity  and  wet-surface  pressure  at  every  time  step.  In  addition,  parallel  files  are 
created  that  retain  restart  information  at  time  intervals  dictated  by  the  user.  The  code 
has  a  variable  time  step  capability  and  can  treat  a  spherical  incident  wave  of  arbitrary 
pressure  profile  and  source  location.  Exponentially  decaying  waves  can  also  be  treated 
by  providing  magnitude  and  decay  information.  In  addition,  incident  wave  pressure  and 
particle  velocity  are  tabulated  and  displayed  with  a  "printer-plot"  package.  If  the 
body  is  in  the  vicinity  of  a  free  surface,  unloading  due  to  reflection  of  the  incident 
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wave  from  the  surface  is  included  and  the  effects  of  bulk  cavitation  on  the  free  field 
pressure  history  are  approximately  treated.  Finally,  selected  response  histories  can 
be  listed  and  then  displayed  for  immediate  examination  using  the "printer-pl ot"  graphics 
package  embedded  both  in  TIMINT  and  in  POSTPR  (see  Sect.  3.8). 

The  computational  strategy  for  the  staggered  solution  procedure  is  embodied  in  the 
following  eight  steps,  assuming  the  solution  is  known  at  time  t: 

(1)  Estimate  the  unknown  structural  restoring  force  K^x  at  t  +  At  from  the  extrapo¬ 
lation  of  current  and  past  values 

(2)  Transform  this  extrapolation  into  fluid  node  values  and  form  the  right-hand 
side  of  the  fluid  equation,  which  also  involves  the  known  incident  pressure 
at  t  +  A  t 

(3)  Solve  the  fluid  equation  and  obtain  a  preliminary  estimate  of  the  total 
pressure  vector  at  t  +  At 

(4)  Transform  fluid  pressures  into  structural  nodal  forces 

(5)  Solve  the  structural  equation  for  the  displacement  and  velocity  vectors 
at  t  +  At 

(6)  Transform  the  computed  structural  restoring  forces  at  t  +  A t  into  fluid 
node  values  and  reform  the  right  hand  side  of  the  fluid  equation 

(7)  Re-solve  the  fluid  equation  and  obtain  refined  values  for  the  total 
pressures  at  t  +  A  t 

(8)  Save  system  responses 

Steps  1,  3,  and  5  constitute  the  basic  staggered  solution  technique,  while  Steps  2  and  4 
are  required  because  of  the  difference  between  the  fluid  and  structural  surface  meshes. 
The  iteration  on  the  fluid  solution  reflected  in  Steps  6  and  7  has  been  added  to  enhance 
accuracy.  Inasmuch  as  the  computation  time  is  overwhelmed  by  the  structural  solution  re¬ 
quirements,  this  requires  only  a  small  increase  in  total  run  time.  The  use  of  a  three- 
point  extrapolation  method  in  Step  1  also  improves  accuracy,  as  discussed  in  [9], 

Implicit  integration  algorithms  have  been  used  for  both  the  fluid  and  structural 
equations.  The  former  is  treated  with  the  3-step  Park  method  [27]  while  the  latter  is 
treated  with  the  "J0M  implementation  of  the  trapezoidal  rule  [28]. 

Typical  input  to  this  processor  includes 


•  Incident  wave  characteristics 


Location  of  source 
Location  of  standoff 
Pressure  profile 
Linear  interpolation 
Cubic  spline  fit 
Exponentially  decaying  wave 

•  Time  step  information 

Start  and  finish  times 
Time  increment  values 

•  Restart  data 

•  Display  directives 

Displacements 

Velocities 

Pressures 

Detailed  user  information  concerning  TIMINT  is  given  in  Appendix  C. 

3.8  THE  RESPONSE  POSTPROCESSOR  POSTPR 

This  utility  is  responsible  for  the  listing  and  "printer-plot"  as  well  as  "vector- 
plot"  graphic  display  of  selected  system  responses  and  pseudo-velocity  shock  spectra. 
Output  files  containing  the  structural  displacement  field  at  user-specified  instants  in 
time  may  also  be  created  from  the  response  history  files  to  provide  "snapshots"  of  the 
deformed  structure.  Some  of  the  same  capabilities  are  also  embedded  in  the  TIMINT  pro¬ 
cessor  for  immediate  selective  scanning  of  the  output.  POSTPR,  however,  is  used  for 
more  detailed  examination  of  the  results  at  a  later  time.  As  a  complete  display  of  all 
structural  and  fluid  DOF  histories  for  even  a  moderate  size  problem  could  run  into 
thousands  of  pages  of  output,  care  must  be  exercised  in  the  selection  of  data  to  be  dis¬ 
played.  Usage  of  this  code  is  discussed  in  Appendix  D. 
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SECTION  IV 
EXAMPLE  PROBLEMS 

This  section  presents  results  generated  by  the  USA  Code  for  three  idealized  under¬ 
water  and  free  surface  shock  problems.  The  structure  studied  in  the  first  problem  is  a 
hollow  circular  beam  of  finite  length,  while  that  involved  in  the  second  and  third  prob¬ 
lems  is  an  infinite,  circular  cylindrical  shell.  In  all  cases,  the  structure  is  excited 
by  a  transverse,  plane  step-wave  of  unit  incident  pressure  and  material  properties  are 
used  that  correspond  to  a  steel  shell  immersed  in  water.  The  input  data  are  normalized 
so  that  the  density  and  speed  of  sound  for  the  fluid  both  equal  unity;  hence,  the  dens¬ 
ity,  Young’s  modulus,  and  Poisson’s  ratio  for  the  structural  material  are  taken  as  7.85, 
98.125,  and  0.3,  respectively.  The  radius  and  wall  thickness  of  the  beam  and  the  cylin¬ 
der  are  1  and  0.01,  respectively,  while  the  length  of  the  beam  is  9.  In  order  to  assess 
the  accuracy  of  the  computational  results,  selected  response  histories  are  compared  with 
those  obtained  by  other  methods. 

4.1  CIRCULAR  BEAM 

The  response  variable  of  primary  interest  in  this  problem  is  the  late-time  asymp¬ 
totic  translational  velocity  of  the  structure.  An  analytical  expression  for  this 
quantity  may  be  obtained  from  (2.5)  by  taking  x  =  j£pg  v(t),  where  is  the  vector  of 
direction  cosines  relating  the  translational  motions  of  the  structural  nodes  and  the 
direction  of  propagation  of  the  plane  incident  wave.  (The  elements  of  Jfps  that  pertain 

to  the  rotational  DOF  are,  of  course,  zero.)  The  introduction  of  this  relation  into  the 

T 

first  of  2.5,  followed  by  premultiplication  of  the  resulting  equation  by  ,  then  yields 

V  =  -is££f<%+£s)  (4>1) 

where  ms  =  Yps  Mflps;  this  follows  from  the  fact  that  =  Kglps  = 

After  the  wave  front  of  the  plane  step-wave  has  enveloped  the  structure,  i.e.,  for 

c  >  V 

£j  =  P  c  uj  I 

£l  =  P  c  UjCtl  -  tA)  (4.2) 

ux  = 

where  U^.  is  the  fluid  particle  velocity  characterizing  the  step-wave,  1_  is  the  unity  vec¬ 
tor,  the  asterisk  denotes  the  temporal  integral  of  the  quantity  beneath  it,  is  the  vec¬ 
tor  of  incident-wave  arrival  times  for  the  fluid  surface  elements,  and  ]rp  is  the  vector 
of  direction  cosines  relating  the  normals  of  the  fluid  elements  to  propagation  vector  of 
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the  plane  incident  wave.  In  addition,  |£g|  <<  |j>g|  for  late-time  motions  (see  Section 
2.2),  so  that  the  second  of  (2.5)  becomes 


-Ps 


V 


t»t 

e 


(4.3) 


The  introduction  of  this  relation  into  (4.1)  then  yields 


(ms  +  ma)  v  =  -  G  +  y^G  M^,  t»t(; 


(4.4) 


T  T  T 

where  the  added  mass  m  =  y^G  M_G  y_  .  But,  from  (4.4),  G  V-,  =  Y„,  so  that  m  is  also 

T  a  -Ps^  ~f~  —Ps  *  ~  — rs  ~P’  a 

given  as  ma  =  YpMfYp- 

T  T  T 

With  £  xps  =  Yp?  the  first  of  (4.2)  yields  Xpg£l  =  p  c  U^y^A^l  -  0,  inasmuch  as 

the  wet  surface  of  the  structure  is  closed.  Hence,  the  right  side  of  (4.4)  vanishes  for 

t  >  t^,  which  gives  the  expected  result  v  =  0„  This  prompts  the  use  of  integrated  forms 

of  (4.1)  and  (4.3)  (with  quiescent  initial  conditions),  which  yields,  instead  of  (4.4), 

* 


(ms  +  ma)  v  =  •  +  t>>t:e 


(4.5) 


The  introduction  of  the  second  and  third  of  (4.2)  into  this  equation  then  provides  the 
desired  expression  for  late-time  asymptotic  translational  velocity 


V 

cx> 


m,  +  m 
d  a 

m  +  m 
s  a 


(4.6) 


where  the  structured  displaced  mass  m^  may  be  shown  to  be  expressible  as  m^  =  p  c  Xp 
A^(tl^  -  Jt^).  Note  that  (4.6)  is  a  general  result,  applicable  to  any  wet-surface  geometry. 


Two  different  uniform  mesh  geometries  were  used  to  study  the  circular  beam.  Ten- 

and  twenty-node  models  were  constructed  with  beam  elements  provided  by  the  structural 

analyzer  SPAR  [29].  The  corresponding  fluid  models  contained  11  and  21  elements  of  equal 

size,  with  12  and  24  circumferential  integration  points  (see  [10]).  In  each  case  there 

was  a  fluid  element  on  each  end  to  account  for  axial  motion  while  all  the  rest  were  evenly 

distributed  along  the  length.  For  the  beam  considered,  m  =  4.439  and  m.  =  28.274;  with 

T  s  d 

m^  determined  as  j  M^Xj  mesh  geometry  has  a  small  effect  on  the  value  calculated  for 

V  .  It  was  found  that  m  =  24.509  for  the  coarse  mesh  and  m  =  24.332  for  the  fine  mesh, 

00  a  a 

which  yield  -  1.823  and  =  1,828,  respectively. 


In  the  response  calculations,  a  constant  time  step  of  0.1  (20  steps  per  envelopment 
period)  was  used  for  both  models;  the  results  are  shown  in  Figures  4-1  and  4-2.  Veloci¬ 
ties  at  the  ends  of  the  beam  are  higher  than  those  at  the  center  because  the  three-dimen¬ 
sional  flow  field  at  the  ends  offers  less  resistance  to  the  plane  wave  excitation  than 
the  two-dimensional  flow  field  at  the  center.  It  is  noted  that  the  responses  of  both 
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TRANSVERSE  VELOCITY  (u) 


TRANSVERSE  VELOCITY  (u) 


models  are  similar  although  those  for  the  finer  mesh  appear  to  tend  to  the  rigid  body 
asymptotic  velocity  more  precisely. 


The  USA(3)  Code  can  also  treat  the  bar  response  of  the  circular  beam  and  to  illus¬ 
trate  this  capability  the  axial  displacement  response  of  the  end  under  the  side-on  load¬ 
ing  is  shown  in  Figure  4-3  for  the  coarse  mesh  problem  only.  The  expected  static  dis¬ 
placement  is  easily  determined  from  the  physical  characteristics  of  the  problem  to  be 
2.293  and  it  is  noted  that  the  computational  result  is  settling  down  to  that  value. 
However,  the  time  required  for  this  to  happen  is  much  greater  than  the  time  for  the  beam 
to  achieve  the  expected  transverse  velocity  shown  in  Figure  4-1.  Thi6  is  because  the 
axial  frequency  of  the  beam  is  low  due  to  the  entrained  fluid.  The  natural  frequency  of 
this  system  can  easily  be  found  by  Rayleigh* s  method  as  follows. 


Under  a  harmonic  excitation  u  =»  B  x  sinw  t  the  maximum  potential  energy  umax  and 

the  maximum  kinetic  energy  T  can  be  shown  to  be 

max 


U 

max 


T 

max 


1  2 
|  B  AEA 

1  2A  32  1  c2  02  2 
24-BApjiu)  +-gBm^£o) 


(4.7) 


where  u  is  the  axial  displacement  measured  from  one  end,  x  is  the  axial  coordinate, 
w  is  the  circular  frequency,  t  is  time,  A  is  the  cross-sectional  area,  E  is  Young* 
modulus,  i  is  the  bar /beam  length  and  mf  is  the  fluid  added  mass  for  the  axial  breath 


ing  mode  of  the  bar /beam.  Equating  U 


max 


to 


T  then  gives  a  value  for  the  period 
max  *  * 


wet 


as 


rwet 


p  ,  /  1  +  3  mT/m 

*ary  f  s 


where 

Pdry  ■  iri.  ✓  p/3E  (4.8) 

m  *  pAi 
s  r 

For  the  problem  at  hand  P^Ty  *  4.617,  *  4.913,  and  mg  *  4.493  which  gives  Pwet  = 

9.597,  about  twice  as  long  as  the  dry  period.  It  is  noted  from  the  peak-to-peak  times 
shown  in  Figure  4-3  that  the  computational  result  is  in  excellent  agreement. 
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AXIAL  DISPLACEMENT  AT  END 


4.0 


TIME  (t) 


Figure  4-3  Axial  Displacement  at  Knd  of  Finite  Beam,  Coarse  Mesh 
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4.2  SUBMERGED  INFINITE  CYLINDRICAL  SHELL 


For  this  problem,  a  72-node,  36-element  SPAR  model  with  a  uniform  circumferential 
mesh  was  constructed.  The  length  of  the  cylindrical  shell  equalled  the  circumferential 
dimension  of  the  square  plate  elements  used  for  the  model;  hence  the  shell  was  one  ele¬ 
ment  long.  Kinematic  constraints  of  zero  axial  displacement  and  no  end  rotation  were 
enforced  through  the  skyline  utility,  as  described  in  Section  3.3.  The  fluid  model  con¬ 
sists  of  36  equally-spaced  elements  around  the  circumference;  the  two-dimensional  nature 
of  the  infinite  shell  geometry  was  simulated  by  exercising  an  option  in  the  fluid  pre¬ 
processor  FLUMAS  that  adds  fictitious  elements  in  the  axial  direction. 

Two-dimensional  n=0,  1,  and  2  modal  response  results  were  generated  by  circum¬ 
ferential  Fourier  decomposition  of  USA-generated  responses.  For  comparison 
DAA  analytical  solutions  were  generated  by  the  method  described  in  [8]  and  [30].  The 
primary  response  variables  of  interest  were  radial  displacement  for  n=0,  radial  and  tan¬ 
gential  velocity  for  n=l,  and  radial  and  tangential  displacement  for  n=2.  A  time  step 
or  0.025  was  used  up  to  t=l;  for  t  between  1  and  2  this  was  increased  0.05,  and  for  t 
greater  than  2  a  time  step  of  0.1  was  used. 

The  USA  and  corresponding  analytical  results  are  shown,  harmonic  by  harmonic,  in  Fig¬ 
ures  4-4  through  4-8.  In  all  cases  the  maximum  errors  made  by  USA  fall  into  the  range  of 
1  to  2%. 

4.3  SEMI-SUBMERGED  INFINITE  CYLINDRICAL  SHELL 

Analytical  solutions  for  the  problem  of  a  shock-wave-excited,  infinite,  circular 
cylindrical  shell  whose  axis  lies  in  the  plane  of  the  free  surface  of  a  semi-infinite 
acoustic  medium  may  readily  be  obtained.  The  appropriate  geometry  for  an  imaging  repre¬ 
sentation  of  this  problem  is  shown  in  Figure  4-9.  An  earlier  treatment  is  described  in 
[31],  where  the  stiffness  of  the  shell  is  neglected  and  an  approximate  expression  is 
obtained  for  the  asymptotic  value  for  vertical  rigid-body  translational  velocity.  In 
this  subsection  it  is  shown  that  both  exact  and  approximate  solutions  for  the  infinite- 
fluid  case  may  be  used  directly  to  obtain  solutions  for  the  case  of  a  semi-infinite  fluid, 
as  follows. 

The  shell  and  fluid  variables  are  first  expanded  in  circumferential  harmonics  as 


'  / 


4-7 


Figure  A-  5  n=l  Radial  Velocity  of  Infinite  Cylinder 
in  Infinite  Fluid 
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Figure  4-8  n*2  Tangential  Displacement  of  Infinite  Cylinder 

in  Infinite  Fluid 
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where  w  and  v  are  the  radial  and  circumferential  shell  displacements,  respectively,  p 
is  the  total  fluid  pressure  and  p^  and  u^.  are  the  pressure  and  radial  fluid-particle 
velocity  of  the  incident  wave.  Because  this  is  a  linear  problem,  the  solution  con¬ 
sists  of  the  superposition  of  two  infinite-fluid  solutions,  the  time-dependent  part  of 
which  will  be  denoted  by  a  superscript  I.  The  result  is 
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(4.10) 


wn(t)  cos  n(0+y) 
(-l)n  iv^(t)  sin  n ( 9+y) 

pJU.t)  cos  n(04-y) 


The  multipliers  (-1)  are  required  if  the  results  of  [3,4,30,32,33]  are  used  to  provide 
the  inf inite-f luid  solutions  for  step-wave  excitation.  This  is  because  of  the  conven¬ 
tion  adopted  there  that  the  incident  wave  first  contacts  the  cylinder  at  e  =  it,  It 
should  be  noted  that  (4.10)  holds  for  any  plane-wave  input  if  the  infinite-fluid  solu¬ 
tion  can  be  found  and  no  cavitation  occurs  in  the  fluid. 

The  trigonometric  terms  in  (4 . 10)can  easily  be  simplified  to  give 


(w(e,t)  | 

^  1 

^  w*(t)  sin  n0 

(v(0,t)  >  = 

2^  (-1)°  sin  ny  ' 

-vl(t)  cos  nG 

|p(r,e,t)J 

n=l 

pT(r,t)  sin  nG 
v  n 
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Note  from  (4. 10)  that  the  solution  for  n=0  is  identically  zero.  Also  note  that  if  the 
time-dependent  part  of  the  free  surface  solution  is  denoted  by  a  superscript  F,  then 
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This  extremely  simple  result  applies  to  exact  solutions  produced  by  the  residual  poten¬ 
tial  method  [3,30,32],  approximate  solutions  produced  by  the  DAA  [3,4],  and  the  approx¬ 
imate  post-envelopment  solutions  given  in  [33]. 

For  a  USA-Code  treatment  of  this  problem,  a  40-node,  19-element  SPAR  half-model  was 
constructed  consisting  of  one  5°  plate  element,  seventeen  10°  plate-elements  and  a 
second  3°  plate-element.  The  angle  y  was  taken  as  45°.  In  other  respects,  this  analy¬ 
sis  paralleled  that  of  the  totally  submerged  cylindrical  shell  described  in  Subsection 
4.2. 

Velocity  responses  for  n=l  are  shown  in  Figures  4-10  and  4-11,  while  displacement  re¬ 
sults  for  n=2  are  shown  in  Figures  4-12  and  4-13.  The  level  of  USA-Code  error  for  n=l 
is  about  1%;  it  is  of  interest  to  note  that  the  approximate  asymptotic  velocity  provided 
by  the  analysis  of  [31]  is  =  2.623,  which  is  about  10%  high. 

The  n=2  post-envelopment  solution  during  the  final  phase  of  shell  envelopment  is 
quite  sensitive  to  temporal  and  spatial  discretization  details.  This  sensitivity  is 
reflected  in  the  level  of  USA-Code  error  exhibited  in  Figures  4-12  and  4-13 ,  and  is 
associated  with  the  pressure  discontinuity  at  the  front  of  the  step-wave.  The  error 
here  is  somewhat  greater  than  that  in  Figures  4-7  and  4-8  ,  as  both  the  incident  wave 
and  its  negative  image  are  now  involved. 

4.4  DAA2  STUDIES 

As  an  initial  test  of  USA-DAA^,  the  infinite-cylinder  case  has  also  been  run 
and  physical  response  variables  have  been  Fourier-decomposed  to  produce  modal 
responses.  For  comparison,  corresponding  exact  solutions  have  been  obtained  analy¬ 
tically  by  the  residual  potential  (RP)  method  [30].  Comparisons  are  shown  in 
Figures  4-14  through  4-18. 

The  difference  between  the  family  of  doubly  asymptotic  solutions  and  the  exact 
solution  shown  in  Figure  4-14  is  intrinsic  to  the  infinite  cylinder  and  occurs 
because  of  the  infinite  value  of  the  n=0  fluid  mass  coefficient.  Hence  the  inter¬ 
action  consists  only  of  the  high  frequency  asymptote  of  the  DAA  family.  For  a 
finite  cylinder,  the  fluid  mass  for  the  breathing  mode  is  always  finite,  so  that 
the  low-frequency  contribution  of  the  DAA  family  does  participate. 
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Velocity  responses  for  n=l  are  shown  in  Figures  4-15  and  4-16,  where  it  is 
noted  that  rr=0  (DAAp  and  n  -  1  provide  rather  wide  bounds  on  the  exact  solution 
in  the  intermediate  time  range,  while  a  value  of  n  *  1/2  is  remarkably  close  to  the 
exact  solution.  In  Figures  4-17  and  4-18,  the  displacement  responses  do  not  show 
such  wide  differences  between  n  -  1/2  and  n  “  1;  however,  they  both  demonstrate  a 
striking  superiority  over  the  DAA^  result.  It  would  appear  from  this  limited 
evidence  that  n  =  1/2  is  probably  optimum  for  the  infinite  cylinder  case. 

In  closing,  it  should  be  emphasized  that  a  study  of  the  stability  of  (2.34) 
and  (2.35)  has  yet  to  be  undertaken.  However,  the  results  obtained  thus  far  are 
very  encouraging  in  this  regard.  No  indication  of  instability  whatsoever  was 
encountered  in  generating  the  infinite  cylinder  results  in  which  the  time  step 
was  doubled  twice  during  the  run,  as  was  done  in  the  DAA^  runs. 
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Figure  4-11  n=l  Tangential  Velocity  of  Infinite  Cylinder  at 
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APPENDIX  A 

USER  INFORMATION  FOR  THE  FLUID  PREPROCESSOR  FLUMAS 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  in  Section  4. 
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THIS  FUNCTIONAL  COMPONENT  OF  THE  UNDERWATER  SHOCK  ANALYSIS  CODE 
CONSTRUCTS  THE  FLUID  MASS  MATRIX  FOR  A  STRUCTURE  SUBMERGED  IN  AN 
INFINITE.  I N  V  I  SC  I D .  INCOMPRESSIBLE  FLJID  B*  THE  BOUNDARY  ELEMENT 
TECHNIQUE.  IT  ;•  ISO  GENERATES  FLUID  MESH  DATA  AND  A  SET  OF 
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QUADRILATERAL 


/ 


The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

The  first  item  needing  explanation  is  the  block  subdivision tabl e .  During  construc¬ 
tion  of  the  mass  matrix  the  code  must  determine  whether  a  particular  fluid  DOF  pertains 
to  a  GEN  (includes  both  GEN  and  CYL  elements)  or  SOR  element.  In  the  latter  case,  it 
must  also  store  the  branch  or  axis  of  the  element,  its  harmonic,  and  also  whether  that 
DOF  corresponds  to  a  cosine  or  sine  function.  DOF  with  similar  characteristics  are 
naturally  kept  together  in  the  same  block.  When  the  mass  matrix  is  automatically  pro¬ 
cessed  in  an  out  of  core  mode  GEN  elements  are  also  partitioned  into  blocks  for  computa¬ 
tional  convenience. 

The  parameters  appearing  in  the  block  subdivision  table  are: 

ISUB  -  block  number 

ITYP  -  GEN  or  SOR 

IBEG  -  first  row  of  block 

IROW  -  number  of  rows  in  block 

IBRA  -  SOR  branch  or  axis 

IHAR  -  harmonic  number 

I FUN  -  COS  or  SIN 

Next,  the  terms  appearing  under  "Fluid  Mesh  Geometric  Arrays"  are  defined  as: 

NCOR  -  number  of  corner  points  for  a  particular  fluid  element 

X,Y,Z  -  global  cartesian  coordinates  of  the  fluid  element  centroidal 

control  point 

NX,NY,NZ  -  components  of  the  outward  unit  normal  vector  for  the  fluid 

element 

NTRA  -  number  of  structural  node  points  that  are  coupled  to  a  par¬ 

ticular  fluid  element  for  the  purpose  of  force  application 

A00, A20, All, A02  -  area  and  moments  and  product  of  inertia  of  fluid  element. 

Used  internally  for  construction  of  the  fluid  mass  matrix 
and  of  the  fluid-structure  transformation  coefficients  for 
general  elements.  For  SOR  elements,  these  values  are  for 
the  sub-elements. 

BII,CII  -  diagonal  terms  of  B  and  C  matrices  used  for  the  construction 

of  fluid  mass  matrix  (see  [14]) 

When  SOR  elements  are  included  in  the  fluid  mesh  the  following  new  terms  will  appear 
in  the  output: 


NSOR  -  number  of  SOR  element 

NFLU  -  DOF  in  fluid  mass  matrix 

RAD  -  radius  of  fluid  element  control  point  from  axis  of  revolution 

NCIR  -  number  of  integration  points  or  sub-elements  used  in  circumferential 
direction 

Local  Fluid-Structure  Transformation  Coefficients  appear  next.  This  is  a  summary 
that  indicates  which  structural  nodes  couple  with  a  particular  fluid  control  point  and 
the  weighting  factor  for  each.  The  weighting  factors  must  always  sum  to  unity  for  any 
fluid  control  point. 

The  geneialized  areas  that  follow  are  simply  A00  for  GEN  elements.  For  SOR  elements 
with  IHAR  =  0  they  are  AQ0*NCIR;  for  all  other  SOR  elements  they  become  . 5*A00*NCIR. 

The  eigenvalues  and  eigenvectors  that  follow  the  listing  of  the  added  mass  matrix 
correspond  to  the  ’’Fluid  Boundary  Mode"  problem  [14].  For  the  infinite  cylindrical 
shell  problem  presented  here,  the  exact  eigenvalues  should  behave  as  1/n  with  correspond¬ 
ing  modes  cos  nQ  and  sin  n0 .  The  first  eigenvalue  listed,  0.11831+04,  is  an  approxima¬ 
tion  to  00  for  n=0  and  it  can  be  seen  that  the  subsequent  eigenvalues  are  relatively  well 
behaved . 

If  a  table  labeled  ’’SUMMARY  OF  1-0  ACTIVITY”  appears  in  the  output,  this  indicates 
that  automatic  out-of-core  processing  has  taken  place.  In  such  a  case  the  ’’Fluid  Boundary 
Mode"  problem  is  not  solved  and  its  diagnostic  characteristics  are  unavailable  to  the  user. 
If  there  are  any  serious  errors  in  the  fluid  mesh  geometry  that  have  remained  undetec*  ad 
through  the  generation  of  the  mass  matrix  these  may  show  up  in  the  construction  of  the 
matrix  [see  Eq.  (2.6)],  i.e.,  the  occurrence  of  factorization  errors  for  the  elements 

in  question. 
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APPENDIX 

USER  INFORMATION  FOR  THE  AUGMENTED 


B 

MATRIX  PREPROCESSOR  AUGMAT 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  in  Section  4. 
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THIS  FUNCTIONAL  COMPONENT  OF  THE  UNDERWATER  SHOCK  ANALYSIS  CODE 
ACCEPTS  DATA  FROM  THE  FLUID  MASS  MATRIX  PROCESSOR  AND  THE 
STRUCTURAL  ANALYZER  TO  CONSTRUCT  THE  SPECIFIC  CONSTANTS  AND  ARRAYS 
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SYMCON  L  TRUE  IF  TRANSLATIONAL  CCNSTRA I  NTS  MUST  BE 

APPLIED  TO  STRUCTURAL  NODES  DUE  TO 
SYMMETRY  CONDITIONS  IF  HALF  OR  QUARTER 
MODELS  ARE  BEING  USED.  OTHERWISE  FALSE. 
CONSTRAINTS  ON  RO  TATI QUA  l  STRUCTURAL 


SHOWN  THAT  THIS  SCALAR  PARAMETER  DOES  HAVE 
A  RELATIONSHIP  WITH  THE  DIAGONAL  LOCAL 
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The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

After  a  summary  of  the  fluid  mesh  geometry  arrays  (see  Appendix  A)  the  first  item 
needing  explanation  is  that  entitled  MFluid  Element  Wetted  Freedom  Indicator'*.  This  is 
simply  a  listing  of  the  input  variable  NDICOS  (see  user  manual)  for  each  fluid  element. 

The  section  "Structural  Grid  Point  Numbers  Associated  With  Internal  Sequence  Num¬ 
bers"  contains  a  correspondence  table  that  relates  the  internal  sequence  numbers  assigned 
by  the  fluid  mass  processor  with  the  external  structural  node  number  assigned  by  the 
user . 

The  next  item  entitled  "Grid  Point  and  Freedom  Number  for  Each  Row  of  Stiffness 
Matrix"  identifies  an  integer  vector  that  is  constructed  by  the  user  in  the  Skyline 
Utility  (see  Figure  3-1,  also  Appendix  E).  For  each  structural  equation  the  entry  in 
the  vector  consists  of  ten  times  the  structural  node  number  plus  the  local  degree  of 
freedom  number. 

The  last  item  requiring  explanation  is  the  "Freedom/Equation  Correspondence  Table". 
This  is  an  integer  matrix  of  6  rows  and  as  many  columns  as  there  are  structural  node 
points.  Any  particular  row  corresponds  to  a  local  degree  of  freedom  number  while  a  col¬ 
umn  corresponds  to  the  internal  sequence  number  for  a  particular  external  node  number. 

The  matrix  entry  for  any  particular  set  of  row  and  column  is  the  structural  equation 
number  for  that  pair. 

Depending  upon  user  input  the  structural  stiffness  matrix  (identifier  "STIF")  may 
then  be  displayed  as  well  as  the  appropriate  fluid  matrices.  The  matrix  called  TMIT 
corresponds  to  D  [see  Eq.  (2.6)],  while  DFDS  denotes  the  sum  of  D  and  Dri  [see  Eq. 
(2.6)].  In  DAA2  runs  Dfl  is  labeled  DAA1  while  Df2  is  labeled  DAA2 . 
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APPENDIX  C 

USER  INFORMATION  FOR  THE  TIME  INTEGRATION  PROCESSOR  TIMINT 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinite  cylindrical  shell  problem  presented  in  Section  4. 
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519  NPREHS  I  NUMBER  OF  FLUID  PRESSURE  HISTORIES  TO  BE 

520  DISPLAYED 

521 

522  NEQHPR  I  FLUID  CONTROL  POINT  NUMBER  FOR  WHICH  A 


The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

First,  the  amount  of  storage  required  for  the  run  given  in  the  output  refers  solely 
to  the  blank  common  that  is  set  in  the  main  program,  UNWASH.  An  error  exit  is  taken  if 
insufficient  storage  is  available  and  the  user  must  see  that  more  is  provided  either  by 
a  recompilation  on  UNIVAC  1100-OS  or  by  a  field  length  request  on  CDC . 

Sector  address  information  for  the  response  and  restart  files  that  is  listed  at 
various  places  in  the  output  is  extremely  important  for  subsequent  restart  runs. 

The  next  item  needing  discussion  is  the  transient  response  tabular  listings.  The 
desired  responses  are  displayed  in  matrix  form  so  that  each  row  contains  the  entire 
history  of  a  particular  degree  of  freedom  except  for  the  first  row  which  is  time.  Each 
column  therefore  contains  the  instantaneous  values  of  the  complete  set  of  response  var¬ 
iables  desired  at  a  particular  time.  Each  row  is  identified  by  the  structural  or  fluid 
node  and  its  degree  of  freedom.  The  letters  D,  V,  and  P  stand  for  displacement,  velocity 
and  pressure,  respectively. 

Although  printer  plots  of  the  transient  response  results  can  be  displayed  as  part 
of  the  run  such  output  has  been  deferred  to  the  post-processing  phase  in  Appendix  D  for 
this  sample  problem. 
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APPENDIX  D 

USER  INFORMATION  FOR  THE  POSTPROCESSOR  POSTPR 


This  appendix  includes  a  copy  of  the  users  manual,  and  a  sample  input  deck  and 
subsequent  output  for  the  infinte  cylindrical  shell  problem  presented  in  Section  4. 
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FORWRT  l  TRUE  IF  PERMANENT  FILES  DENOTED  BY  XVPNAM 


WERE  CREATED  USING  UNFORMATTED  FORTRAN 
WRITE.  OTHERWISE  FILES  WERE  CREATED  BY 
DIRECT  TRANSFER  USING  THE  DATA  MANAGEMENT 

system  dmgasp 
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The  following  discussion  is  provided  as  an  aid  to  user  understanding  of  the  sample 
output  that  is  included  here. 

The  input  deck  shown  on  the  next  page  requests  vector  plots  for  both  the  transient 
response  histories  and  pseudo-velocity  shock  spectra.  This  is  appropriate  if  the  DISSPLA 
plot  package  is  available  at  the  us*  installation.  Otherwise  appropriate  modifications 
must  either  be  made  to  use  a  different  plot  package  or  the  input  deck  should  be  modified. 
In  any  case  the  printer  plot  package  is  resident  in  USA. 

The  format  used  for  listing  the  pseudo-velocity  shock  spectra  is  similar  to  that 
used  for  the  display  of  the  transient  response  histories  shown  in  Appendix  C  except  that 
the  first  row  is  now  frequency  rather  than  time. 
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DISPLACEMENT  RESPONSE  OF  STRUCTURAL  NODE  19,  FREEDOM  NUMBER 
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3.0 
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VELOCITY  RESPONSE  OF  STRUCTURAL  NODE  19,  FREEDOM  NUMBER 


2.04 


VELOCITY  RESPONSE  OF  STRUCTURAL  NODE  37,  FREEDOM  NUMBER 
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PRESSURE  RESPONSE  OF  FLUID  NODE 


ENO  DISSPLA  —  22151  VECTORS  GENERATED -IN  11  PLOT  FRAMES. 
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PSEUDO-VELOCITY  SHOCK  SPECTRUM  FOR  STRUCTURAL  NOOE  1,  FREEDOM  NUMBER 


PSEUDO-VELOCITY  SHOCK  SPECTRUM  TOR  STRUCTURAL  node  19,  FREEDOM  NUMBER 


PSEUDO— VE LOC I T Y  SHOCK  SPECTRUM  FOR  STRUCTURAL  NOQE  37.  FREEDOM  NUMBER 


APPENDIX  E 

USER  INSTRUCTIONS  FOR  INTERFACING  WITH  USA 


To  use  the  Underwater  Shock  Analysis  (USA)  Code  in  its  linear  stand-alone  mode,  the 
user  must  first  construct  a  permanent  data  file  that  contains  the  structural  mass  and 
stiffness  matrices  and  some  assorted  bookkeeping  information.  The  purpose  of  this  appen¬ 
dix  is  to  describe  the  structure  of  the  file  and  to  specify  how  it  is  to  be  created.  At 
this  time  utility  routines  that  carry  out  this  task  have  been  written  for  SPAR,  NASTRAN 
and  GENSAM.  An  abbreviated  form  of  this  file  is  also  required  when  USA  is  coupled  with 
a  non-linear  structural  analyzer  and  such  an  interface  also  exists  for  STAGS. 

USA  contains  the  data  management  utility  module  DMGASP  that  carries  out  all  data 
transfer  activities  between  core  and  peripheral  storage.  This  is  done  by  unformatted 
and  unbuffered  data  transmissions  and  it  is  imperative  that  DMGASP  be  used  to  create  the 
structural  interface  file.  Otherwise  the  user  must  supply  or  have  access  to  a  similar 
means  of  direct  transfer.  Section  3  of  [19]  contains  a  comprehensive  discussion  of  the 
half-dozen  or  so  DMGASP  commands  that  are  required  to  activate,  position,  write  upon, 
read  from,  and  free  a  peripheral  storage  device.  Subsidiary  commands  also  exist  for 
error  handling  and  listing  of  selected  information  pertaining  to  auxiliary  storage. 

The  current  configuration  of  USA  uses  a  diagonal  mass  matrix  associated  with  a 
lumped  mass  representation  of  the  structure,  assumes  that  there  is  no  velocity  dependent 
structural  damping,  and  further,  only,  single  precision  matrices  may  be  processed.  In 
addition,  if  the  stiffness  matrix  has  been  reordered  or  reduced  in  any  way  for  input 
to  USA  the  mass  matrix  must  also  be  reordered  or  condensed  so  that  its  degrees  of 
freedom  (DOF)  are  the  same  and  appear  in  the  same  order  as  in  the  stiffenss  matrix. 
Finally,  the  stiffness  matrix  must  be  placed  in  a  multi-block*  skyline  format  as  dis¬ 
cussed  in  the  SKYPUL  manual  [23].  This  description  consists  of  a  Matrix  Master  Record 
(MMR)  followed  by  a  series  of  Matrix  Value  Records  (MVR)  which  contain  the  numerical 
values  of  the  matrix.  These  are  the  only  constructs  the  user  need  be  concerned  with; 
all  others  required  are  already  embedded  in  USA.  During  construction  of  the  MMR  a  logi¬ 
cal  device  index  (LDI)  must  be  set  in  the  record  which  USA  will  access  later.  For  UNIVAC 
operation,  this  should  be  set  equal  to  twenty  (20),  while  for  CDC  operation  this  should 
be  set  as  two  (2). 

The  USA  code  now  includes  four  routines  in  the  AUGMAT  processor  that  can  be  used  to 
facilitate  assembly  of  the  MMR  and  the  MVR  for  any  symmetric  matrix.  The  routine  BLKSKY 
is  used  to  sequentially  call  SETPAR  which  determines  the  skyline  profile,  SETMMR  which 

*  For  small  problems  a  single  block  is  permissible. 


E-l 


constructs  the  matrix  master  record  and  finally  SETMVR  which  produces  the  matrix  value 
records.  However,  the  user  must  ensure  that  the  matrix  can  be  provided  row  by  row  (or 
column  by  column)  with  any  and  all  zeroes  filled  in.  Please  contact  the  developer  for 
further  details. 

The  file  structure  required  is  shown  in  Table  E-l  where  NDOF  stands  for  the  number 
of  structural  DOF  which  USA  must  process.  NMMR  is  the  number  of  words  in  the  matrix 
master  record,  and  NWBL  is  the  number  of  words  in  each  matrix  value  record  (which  is  ex¬ 
pected  to  be  the  same  for  each  record).  NWBL  should  also  be  an  integer  multiple  of  448 
for  most  efficient  use  of  storage. 

The  Grid  Point/DOF  vector  consists  of  an  integer  value  for  each  global  DOF  from  1 
through  NDOF  that  is  constructed  as  ten  times  the  external  node  number  plus  the  local 
DOF  number  that  apply  to  that  particular  structural  equation. 

For  example,  if  the  87th  DOF  to  appear  in  the  mass  and  stiffness  matrices  corres¬ 
ponds  to  the  second  degree  of  freedom  at  a  node  identified  externally  as  4637  then  the 
87th  entry  in  the  Grid  Point/DOF  vector  would  be  46372.  Local  translational  degrees  of 
freedom  should  be  numbered  1-3,  rotational  degrees  of  freedom  should  be  numbered  from 
4-6  and  any  others  should  be  numbered  with  7-9.  If  more  than  9  degrees  of  freedom  are 
carried  at  any  node  it  is  a  simple  matter  to  change  the  factor  of  ten  to  one  hundred  in 
a  few  places  in  USA  to  accommodate  this. 

It  should  be  noted  that  records  1-4  are  always  accessed  by  the  USA  pre-processor  AUGMAT 
before  the  time  integration  phase  of  the  analysis  commences?  This  portion  of  the  file 
is  required  for  both  USA  in  the  linear  stand-alone  mode,  and  for  USA  when  it  is  inter¬ 
faced  with  a  nonlinear  structural  analyzer.  In  this  latter  case,  the  fifth  and  succeed¬ 
ing  records  do  not  exist. 

There  is  a  minor  difference  in  the  way  that  the  fifth  and  succeeding  records  are 
constructed  in  the  SPAR  and  NASTRAN  utilities.  In  the  SPAR  interface  one  pass  is  made 
through  the  stiffness  matrix  to  first  determine  its  connectivity  and  set  up  the  book¬ 
keeping  to  construct  the  MMR.  The  MMR  is  then  written  and  a  second  pass  is  made  through 
the  matrix  to  write  the  MVRTs  immediately  following  the  MMR.  In  the  NASTRAN  interface 
only  one  pass  is  made  through  the  stiffness  matrix.  Hence,  the  MMR  is  not  constructed 
until  all  the  MVRTs  have  been  written.  To  follow  the  order  required  by  Table  E-l  space 
is  left  in  the  file  for  the  MMR  to  be  written  in  its  proper  order  on  the  file  after  the 
entire  set  of  MVR's  have  been  written.  Because  of  this  there  will  generally  be  a  buffer 
area  of  irrelevant  data  (garbage)  between  the  MMR  and  the  first  MVR  on  a  NASTRAN  file 
in  contrast  to  a  SPAR  file. 

*  The  fifth  and  succeeding  records  can  be  accessed  in  AUGMAT  if  the  user  wishes  to 
check  the  stiffness  matrix. 
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Diagonal  Mass  Matrix 


NDOF 


Grid  Point/DOF  Vector 


NMMR 


Matrix  Master  Record  for 
Stiffness  Matrix 


First  Matrix  Value  Record 
for  Stiffness  Matrix 

Second  Matrix  Value  Record 
for  Stiffness  Matrix 


Last  Matrix  Value  Record 
for  Stiffness  Matrix 


PI 


The  amount  of  space  currently  allowed  for  the  MMR  on  a  NASTRAN  skyline  file  is  10 
PRU's  (640  words)  on  CDG  and  10  sectors  (280  words)  on  UNIVAC.  These  values  translate 
into  a  current  limit  of  283  and  80  skyline  blocks  for  the  two  systems,  respectively. 

If  more  capability  is  desired  the  statement  MMRPRU  -  10  at  the  beginning  of  subroutine 
KDD  of  the  NASTRAN  skyline  utility  can  be  increased  to  suit  the  user’s  needs. 

The  SK.YPUL  processor  has  the  ability  to  apply  constraints  due  to  symmetry  or 
attachment  to  ground  during  the  time  integration.  Structural  DOF  that  are  to  remain 
zero  must  have  their  associated  diagonal  location  pointers  (LDP)  flagged  with  a  negative 
sign  during  construction  of  the  MVR  and  it  is  highly  recommended  that  this  capability 
be  included  in  every  USA  structural  interface  utility. 
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ATTN: 

Code  1740.5 

ATTN: 

Code  1740.6 

Marine  Corps 
Department  of  the  Navy 
ATTN:  POM 

Marine  Corp  Dev  &  Ed  Cmd 
Department  of  the  Navy 

ATTN:  D091,  J.  Hartneady 

Naval  Civil  Engrg  Lab 

ATTN:  Code  L51,  J,  Crawford 

Naval  Coastal  Systems  Lab 
ATTN:  Code  741 
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Naval  Research  Lab 


ATTN: 

Code  8440,  G. 

O'Hara 

ATTN: 

Code  8445 

ATTN: 

Code  8301 

ATTN: 

Code  8403,  R. 

Bel  sham 

ATTN: 

Code  6380 

ATTN: 

Code  8100 

ATTN: 

Code  8406 

ATTN: 

Code  8404,  H. 

Pusey 

ATTN: 

Code  2627 

Sea  Systems  Cmd 

ATTN: 

SEA-3221 

ATTN: 

SEA-09G53 

ATTN: 

SEA-033 

ATTN: 

SEA-06J,  R.  Lane 

ATTN: 

SEA-0351 

ATTN: 

SEA-323 

ATTN: 

SEA-9931G 

ATTN: 

SEA-08 

Naval  Surface  Weapons  Ctr 
ATTN:  Code  R14 
ATTN:  Code  R13 

ATTN:  Code  U401,  M.  Kleinerman 

ATTN:  Code  F34 

ATTN:  Code  R15 

ATTN:  Code  F31 

ATTN:  Code  R14 

ATTN:  Code  RIO 

Naval  War  College 

ATTN:  Code  E-ll,  Tech  Svcs 

Naval  Weapons  Ctr 

ATTN:  Code  266,  C.  Austin 
ATTN:  Code  233 
ATTN:  3263,  J.  Bowen 


Naval  Electronic  Systems  Cmd 
ATTN:  PME  117-21 

Naval  Electronics  Systems  Cmd 
ATTN:  Commander 

Naval  Explosive  Ord  Disposal  Fac 

ATTN:  Code  504,  J.  Petrousky 

Naval  Facilities  Fngrg  Cmd 
ATTN:  Code  04 B 

Naval  Material  Cmd 

ATTN:  MAT  08T-22 

Naval  Ocean  Systems  Ctr 
ATTN:  Code  4471 
ATTN:  Code  013,  E.  Cooper 

Naval  Postgraduate  Sch 
ATTN:  Code  69NE 
ATTN:  Code  1424,  Lib 

Naval  Surface  Weapons  Ctr 
ATTN:  W.  Wishard 
ATTN:  Tech  Lib  &  Info  Svcs  Br 


Naval  Weapons  Evaluation  Fac 
ATTN:  R.  Hughes 
ATTN:  Code  10 
ATTN:  G.  Binns 
ATTN:  Code  210 

Naval  Weapons  Support  Ctr 

ATTN:  Code  70553,  D.  Moore 

Naval  Underwater  Systems  Ctr 

ATTN:  Code  401,  J.  Kalinowski 
ATTN:  Code  401,  J.  Patel 

Naval  Underwater  Systems  Ctr 

ATTN:  Code  363,  P.  Paranzino 
ATTN:  Code  EM 

Office  of  Naval  Research 

ATTN:  Code  474,  N.  Per rone 

Strategic  Systems  Proj  Ofc 

Department  of  the  Navy 
ATTN:  NSP-272 
ATTN:  NSP-273 
ATTN:  NSP-43 
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Office  of  the  Chief  of  Naval  Ops 
ATTN:  OP  982 
ATTN:  OP  03EG 
ATTN:  OP  65 

ATTN:  OP  654C3,  R.  Piacesi 

ATTN:  OP  981N1 

ATTN:  OP  981 

ATTN:  OP  21 

ATTN:  OP  957E 

ATTN:  OP  987 

ATTN:  OP  37 

ATTN:  OP  982E,  M.  Lenzini 

ATTN:  OP  951 

ATTN:  OP  225 

ATTN:  OP  60505 

ATTN:  OP  953 

ATTN:  OP  098T8 

ATTN:  OP  223 

DEPARTMENT  OF  THE  AIR  FORCE 


Rome  Air  Development  Ctr,  Air  force  Sys  Cmd 
ATTN:  Commander 
ATTN:  RBES,  R.  Mair 
ATTN:  TSLD 

Strategic  Air  Command 

Department  of  the  Air  force 
ATTN:  NR I ,  STINFO  Lib 

United  States  Air  Force  Academy 
ATTN:  DFCEM,  W.  Fluhr 

DEPARTMENT  OF  ENERGV 

Department  of  Energy 

Albuquerque  Ops  Ofc 
ATTN:  CTID 

Department  of  Energy 
ATTN:  OMA/RD&T 


Air  Force  Institute  of  Technology 
ATTN:  Commander 
ATTN:  Lib 


Department  of  Energy 
Nevada  Ops  Ofc 

ATTN:  Mail  &  Rees  for  Tech  Lib 


Air  Force  Ofc  of  Scientific  Rsch 
ATTN:  NA,  B.  Wolfson 

Air  Force  Systems  Cmd 
ATTN:  DLW 
ATTN:  R.  Cross 

Air  Force  Weapons  Lab 
Air  Force  Systems  Command 
ATTN:  NTED 
ATTN:  SUL 

ATTN:  NTE,  M.  Plamondon 
ATTN:  NTES-C,  R.  Henny 
ATTN:  NTES-G,  S.  Melzer 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
ATTN:  IN 

Ballistic  Missile  Ofc 
Air  Force  Systems  Cmd 
ATTN:  DEB 


OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 
ATTN:  OSR/SE/F 

Department  of  the  Interior 

U.S.  Geological  Survey 
ATTN:  0.  Roddy 

Federal  Emergency  Mgt  Agency 

National  Sec  Ofc  Mitigation  &  Rsch 
ATTN:  Mitigation  &  Rsch  Div 

NASA 

Ames  Research  Ctr 

ATTN:  R.  Jackson 
ATTN:  F.  Nichols 

U.S.  Nuclear  Regulatory  Com 

ATTN:  R.  Whipp  for  Div  Sec  *  L.  Shao 

DEPARTMENT  OF  ENERGY  CONTRACTORS 


Deputy  Chief  of  Staff 
Research,  Development  &  Acq 
Department  of  the  Air  Force 
ATTN:  R.  Steere 
ATTN:  AFRDQI 

Deputy  Chief  of  Staff 
Logistics  &  Engrg 
Department  of  the  Air  Force 
ATTN:  LEEE 


Los  Alamos  National  Lab 
ATTN:  R.  Sanford 
ATTN:  G.  Spillman 
ATTN:  R.  Whitaker 
ATTN:  A.  Davis 
ATTN:  MS  364 
ATTN:  M/S634,  T.  Dowler 
ATTN:  MS  670,  J.  Hopkins 
ATTN:  J.  Kodis 
ATTN:  M.  Sandford 


Foreign  Technology  Oiv  Oak  Ridge  National  Lab 

Air  Force  Systems  Command  ATTN:  Civil  Def  Res  Proj 


ATTN 

TQTD 

ATTN:  Central  Rsch  Lib 

ATTN 

NIIS,  Lib 

ATTN 

SDBG 

Sandia  Labs 

ATTN 

SDBF,  S.  Spring 

Livermore  Lab 

ATTN:  Lib  &  Scty  Class  Div 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Hydro-Conduit  Corp 

ATTN:  J.  Duncan 


Kaman  AviDyne 
ATTN:  i 
ATTN:  I 
ATTN:  1 
ATTN:  I 


G.  Zartarian 
N.  Hobbs 
Lib 

R.  Ruetenik 


Kaman  Sciences  Corp 


lib 

F.  Shelton 
D.  Sachs 


Kaman  Tempo 

ATTN:  DAS  I AC 

Kaman  Tempo 

ATTN:  J.  Moulton 

Karagozian  and  Case 

‘  ATTN:  J.  Karagozian 


Sandia  National  Lab 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN:  L.  Vortman 

ATTN 

Tech  Info  Ctr 

ATTN:  3141 

ATTN 

B.  Almroth 

29  cy  ATTN 

T.  Geers 

Lawrence  Livermore  National  Lab 

4  cy  ATTN 

J.  DeRuntz 

ATTN:  S.  Erickson 

4  cy  ATTN 

C.  Felipta 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN 

TIC,  Lib 

BDM  Corp 

ATTN:  T.  Neighbors 

M  &  T  Company 

ATTN:  Corporate  Lib 

ATTN 

D.  McNaight 

ATTN:  A.  Lavagnino 

Management  Science  Associates 

Bel  1  Telephone  Labs 

ATTN 

K.  Kaplan 

ATTN:  J.  White 

McDonnell 

Douglas  Corp 

California  Institute  of  Technology 

ATTN 

R.  Halprin 

ATTN:  T.  Ahrens 

University  of  Miami 

California  Research  &  Technology,  Inc 

ATTN 

S.  Wang 

ATTN:  M.  Rosenblatt 

ATTN 

F.  Tappert 

ATTN:  Lib 

ATTN 

B.  Lemehaute 

ATTN:  K.  Kreyenhagen 

ATTN:  S.  Schuster 

NKF  Engineering  Associates,  Inc 

ATTN 

R.  Belsheim 

University  of  California  at  San  Diego 

ATTN:  W.  Van  Dorn 

Pacific-Sierra  Rsch  Corp 

ATTN 

H.  Brode 

Columbia  University 

ATTN:  F.  Dimaggio 

Pacifica  Technology 

ATTN:  H.  Bleich 

ATTN 

R.  Bjork 

ATTN 

G.  Kent 

University  of  Denver 

ATTN 

R.  Allen 

ATTN:  U.  Wisotski 

Physics  Appl ications ,  Inc 

Electric  Power  Rsch  Institute 

ATTN 

C.  Vincent 

ATTN:  G.  SI  iter 

Physics  International  Co 

General  Dynamics  Corp 

ATTN 

J.  Thomsen 

ATTN:  M.  Pakstys 

ATTN 

E.  Moo^e 

ATTN:  J.  Mador 

ATTN 

Tech  Lib 

ATTN:  J.  Miller 

ATTN 

F.  Sauer 

ATTN 

L.  Behrmann 

General  Electric  Co 

ATTN:  M.  Bortner 

University  of  Pittsburgh 

ATTN 

M.  Willi ms ,  Jr 

R  &  D  Associates 

ATTN:  P.  Haas 

Science  Appl ications ,  Inc 
ATTN:  Tech  Lib 
ATTN:  H.  Wilson 
ATTN:  M.  Mkcay 

Southwest  Research  Institute 
ATTN:  W.  Baker 
ATTN:  A.  Wenzel 

SRI  International 

ATTN:  G.  Abrahamson 
ATTN:  W.  Wilkinson 
ATTN:  A.  Florence 

Teledyne  Brown  Engineering 
ATTN:  J.  Ravenscraft 
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Systems,  Science  &  Software,  1 

ATTN 

T.  McKinley 

ATTN 

D.  Grine 

ATTN 

T.  Cherry 

ATTN 

Lib 

ATTN 

T.  Riney 

ATTN 

R.  Sedgewick 

ATTN 

K.  Pyatt 

Tetra  Tech 

Inc 

ATTN 

L.  Hwang 

TRW  Defense  &  Space  Sys  Group 

ATTN 

F.  Pieper 

ATTN 

P.  Dai 

ATTN 

G,  Hulcher 

ATTN 

E.  Wong 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Defense  &  Space  Sys  Group 
ATTN:  P.  Bhutta 
ATTN:  W.  Lipner 
ATTN:  A.  Narevsky 
ATTN:  D.  Jortner 
ATTN:  B.  Sussholtz 
ATTN:  A.  Feldman 
ATTN:  Tech  Info  Ctr 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  Isenberg 
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